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Abstract By direct Monte-Carlo simulations it is shown that the anisotropic flows can be 
£NJ ■ successfully described in the colour string picture with fusion and percolation provided anisotropy 

of particle emission from the fused string is taken into account. Quenching of produced particles 
in the strong colour field of the string is the basic mechanism for this anisotropy. The concrete 
realization of this mechanism is borrowed from the QED. Due to dependence of this mechanism on 
the external field strength the found flows grow with energy, with values for V2 at LHC energies 
greater by 15% than at RHIC energies. 

(N ! 1 Introduction 

> ■ 

. The study of anisotropic flows in the spectra of particles produced in high-energy heavy-ion collisions 

proves to be very informative as to the dynamics of the underlying emission mechanisms. Observed 
experimental values for the anisotropic flow coefficients v n impose strong restrictions on the models 
which try to describe the data. Of special importance are the coefficients with odd values of 
n, which disappear on the average and can only be observed on the event-by-event basis. They 
carry information about the details of the evolution of the hot nuclear matter in the overlap and 
its fluctuations during the collision. These fluctuations in principle may come from both initial 
conditions for the nuclear matter evolution and dynamics of the successive evolution itself. Simple 
phenomenological initial conditions based on fluctuations in the distribution of participant nucleons 
followed by the hydrodynamical evolution allow to describe the third harmonic coefficient V3 quite 
successfully [SI El IT0| ITT] . Remarkably less phenomenological initial conditions which follow from 
the preasymptotic evolution of emitted gluons in the Colour Glass Condensate approach seriously 
underestimate V3 [5J |TT] unless this evolution is treated in a rigorous non-perturbative manner and 
the subsequent evolution has a non- vanishing viscosity [12] 

In [13] we drew attention to the fact that the elliptic flow could also be naturally explained 
in the colour string approach, with string fusion and percolation [HI [15j [161 H7], quite successful 
in the description of particle production and correlations in the soft part of spectra [181 [191 120] . 
In this approach partons emitted at some point have to pass a certain length before they appear 
outside and are observed. On their way they cross a strong colour field inside the strings and emit 
gluons, so that their energy diminishes. As a result, the observed particle energy turns out to be 
smaller than at the moment of its creation and this energy quenching depends on the path length 
passed inside the nuclear matter and so on the direction of its transverse momentum p. Very crude 
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estimates made in |13| confirmed that a sizable elliptic flow followed from this mechanism. Its 
centrality and transverse momentum dependence qualitatively agreed with the behaviour of the 
RHIC data [2H E2J [231 IM] • These results also supported the ones obtained in a similar framework 
using different simplified methods [251 ESI [27] . 

It was stressed in [T3] that any trustworthy and comprehensive quantitative results could only 
be achieved in more elaborate studies based on detailed Monte-Carlo simulations. This is especially 
true for odd flow harmonics which, as mentioned, can only be studied on the event-by-event basis, 
with all sources of fluctuations taken into account. Calculations along these lines constitute the 
subject of this paper. On the general they confirm crude predictions made in [T3] and also allow 
for a comprehensive study of azimuthal anisotropy in particle production and thus determination of 
higher anisotropic coefficients v n with n > 2. Note that the colour string approach combines both 
the initial conditions created in the course of collisions and the subsequent evolution of the nuclear 
matter in the collision zone modeled by the fusion of produced strings and their final decay into 
observed hadrons. Fluctuations may be present at all stages of this process. Our calculations show 
that the decisive role for the successful description of both V2 and v% is played by the fluctuations 
in the initial geometry of the collision plus those in the string fusion. Fluctuations in the final 
production of observed particles play a minor role. Our model has a single new parameter, the 
quenching coefficient, which characterizes energy loss of the produced parton in the string matter. 
We adjust this parameter to agree with coefficient V2 observed in mid-central Au+Au collisions 
at 200 GeV and integrated over the transverse momenta. With thus adjusted parameter we are 
able to explain V2 at all centralities and transverse momenta for collisions at 62.4 GeV, 200 GeV 
and 2.76 TeV. The third harmonic V3 comes out a little smaller than the average observed values. 
However the event-by-event fluctuations in ^3 (as well as in V2) are calculated to be quite large so 
that the observed values are found to be within the calculated ones when these fluctuations are 
taken into account. 

The structure of the paper is as follows. In the next section we briefly discuss our model and 
also present our quenching mechanism borrowed from QED. In Section 3 we describe our Monte- 
Carlo procedure to study the anisotropic flow on the event-by-event basis. Section 4. presents our 
numerical results. In Conclusion we draw some lessons from our study and point out some possible 
further refinements. 



2 The model 

One needs to combine two ingredients to have anisotropic flows. First, string fusion has to generate 
clusters which are azimuthally asymmetric and emit particles anisotropically. But by itself this will 
not produce any elliptic flow unless the distribution of these clusters in the transverse plane is also 
azimuthally asymmetric. This latter phenomenon can only occur if the clusters are large enough 
to feel the asymmetric form of the overlapping region in the collision. Such clusters arise in the 
process of percolation of fused strings. 

In our model it is assumed that at the moment of the collision color strings are stretched between 
partons of the colliding nuclei. Since they are many and so overlap in the transverse space, they 
fuse and percolate to form macroscopic clusters at some critical string density p = p c where 



rio is the transverse area of simple strings, N is their number and f2 is the nuclear overlap area. 
Starting from the moment of their formation strings decay into particles (quark-antiquark pairs), 
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which process we describe using the well-known formalism for pair creation in a strong electromag- 
netic field. According to this mechanism, in its simplest version, the particle distribution at the 
moment of its production by the string is 

P(p,4>) = Ce-~. (2) 

where po is the particle initial transverse momentum, T is the string tension (up to an irrelevant 
numerical coeffcient) and C is the normalzation factor. However po is different from the observed 
particle momentum p because the particle has to pass through the fused string area and emit gluons 
on its way out. So in fact in Eq. ([2]) one has to consider po as a function of p and path length 
I inside the nuclear overlap: pq = f(p,l(4>)) where <j) is the azimuthal angle. Note that Eq. ([2]) 
describes the spectra only at very soft pq. To extend its validity to higher momenta one may use 
the idea that the string tension fluctuates, which transforms the Gaussian distribution into the 
thermal one [291130]: 

PQ 

P{p,<j>) = Ce vWi. (3) 

Radiative energy loss has been extensively studied for a parton passing through the nucleus or 
quark-gluon plasma as a result of multiple collisions with the medium scattering centers [31] . In our 
case the situation is somewhat different: the created parton moves in the external gluon field inside 
the string. In the crude approximation this field can be taken as being constant and orthogonal 
to the direction of the parton propagation. In the same spirit as taken for the mechanism of pair 
creation, one may assume that the reaction force due to radiation is similar to the one in the QED 
when a charged particle is moving in the external electromagnetic field. This force causes a loss of 
energy which for an ultra-relativistic particle is proportional to [its momentum x field] 2 / 3 [32j: 

^ = -0.12e 2 (e^)) 2/3 , (4) 
where E is the external electric field. Eq. @ leads to the quenching formula 

Po(p,0 =p(l + Kp- 1 / 3 T 2 / 3 lf, (5) 

where we identified eE/n = T as the string tension. As mentioned in the Introduction, the 
quenching coefficient k is adjusted to give the experimental value for the coefficient V2 in mid- 
central Au+Au collisions at 200 GeV, integrated over the transverse momenta. 

Of course the possibility to use electrodynamic formulas for the chromodynamic case may raise 
certain doubts. However in [33J it was found that at least in the N = 4 SUSY Yang-Mills case 
the loss of energy of a coloured charge moving in the external chromodynamic field was given by 
essentially the same expression as in the QED. 

Note that from the moment of particle creation to the moment of its passage through other 
strings a certain time elapses depending on the distance and particle velocity. During this time 
strings decay and the traveling particle will meet another string partially decayed, with a smaller 
colour Q than at the moment of its formation. So one has to consider a non-static string distribution 
with string colours evolving in time and gradually diminishing until strings disappear altogether. 
The time scale of this evolution is estimated to be considerably greater than time intervals charac- 
teristic for partons traveling inside the string matter. However the effect of string decay with time 
is noticeable and we take it into account in our calculations. 
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To study the time evolution of strings we again turn to the Schwinger mechanism. For it one 
has the probability of pair creation in unit time and unit volume as [28J 

T vt = ±-T*e-i, (6) 

where again T stands for eE/ir in QED. For a realistic string the volume V = SL Z where S is the 
string transverse area and L z is the longitudinal dimension of the string. For the string of colour 
Q is T = QTq where Tq is the string tension of the ordinary string with Q = 1. The average 
transverse momentum squared of the emitted quark-antiquark pair < p\ > is just T.To estimate 
L z we assume that the string emits a pair when its energy is equal to 2 < po >= y/T, which gives 
L z = 1/ y/T, so that we get the average probability in unit time 

r t = — T 3 / 2 S. (7) 
Air 

The string colour diminishes by unity with each pair production. So we find an equation which 
describes the time evolution of the string colour Q(t) 

^ = -«Q 3/2 (t) (8) 

with the solution 

(1 + jiaWQo) 2 

where Qo is the initial colour at the moment of the string creation. Coefficient a = T S/(2tt) 
depends on the string transverse area S. As will be explained in next section, we use the picture in 
which the fused string is in fact modeled by a set of " ministrings" formed at intersections of simple 
strings with the same area as the simple string, but greater color. This gives a = 0.03 1/fm. This 
value has been used in our calculations. The average color of ministrings is of the order 2-3. So 
it changes only by 30 -50 % even when the emitted parton travels 5 fm of distance. Still this effect 
is felt in the calculations. In fact it practically does not change the results but changes the value 
of the quenching coefficient k, which in any case is to be adjusted, as explained above. In this 
sense our results are practically independent of the concrete choice of a in the reasonable interval 
of values. 



3 Monte-Carlo simulations 

3.1 Generalities 

In principle Monte-Carlo simulations of the quenching in the fusing string scenario seem to be 
straightforward. One models strings in the nuclear overlap by discs of a given radius r^. In an event 
N discs are assumed to be distributed in the transverse plane in agreement with the distribution of 
the participant ("wounded") nucleons, given by product of the nuclear profile functions. Inside the 
transverse area of each pair of colliding nucleons strings may be taken to be distributed according 
to the nucleon density of the Gaussian form. The number of discs is to be chosen in agreement 
with the observed value of the percolation parameter p given by Eq. ([1]) Values of the percolation 
parameter p can be taken from [M] for Au-Au collisions at 62.4 and 200 GeV. They were extracted 
from the observed distributions in the transverse momentum at different centralities and so depend 
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Figure 1: Percolation parameter as a function of impact parameter for Au-Au collisions at 62.4 
(lower curve) and 200 GeV borrowed from [M] 



on the impact parameter b. They are shown in Fig. [TJ For the LHC energy E = 2.76TeV we used 
the conclusions in [35 1 136] that the values of p are roughly 4 times larger than at RHIC energies. 

The number of particpant nucleons for Au-Au collisons at 200 GeV and given centralities, or 
equivalently, b can be borrowed from [11] , They are reproduced in Table 1. 

Taking the radius of both colliding nuclei R = A 1 ^Rq with Rq = 1.2 fm and ro = 0.3 fm, for 
Au — Au collisions at 200 GeV one gets from Eq. ([I]) and Fig. [T] the number of strings given at 
a given b. Comparison with Table 1. then shows that for Au-Au collisions at 200 GeV each pair 
of participant nucleons gives rise to approximately 7 colour strings. Assuming that the number 
of participant nucleons is purely geometrical and does not change significantly with energy we 
conclude that for Au+Au collisions at 62.4 GeV and 2.76 Tev the numbers of strings per each pair 
of participant nucleons are 4 and 28 respectively. 

The strings modeled by discs may overlap in the transverse area and form clusters of different 
number n of fused strings and form. Observed particles are emitted from each cluster with the 
average multiplicity 



/ nVt 



and average transverse momentum squared 
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(10) 



(11) 



for the k-th cluster of n fused strings. Here £l n k is the transverse area of the cluster and an d 
Pq are the multiplicity and transverse momentum squared for a simple string. In our picture each 
particle emitted from a given point in a cluster has to pass a certain path in the overlap area 
before being observed. A part of it has to pass through the same or different clusters and so looses 
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Table 1: Number of particpant nucleons for given centrality and b for Au+Au at 200 GeV 



Centrality 


6 (fm) 


Npart 


0-5 % 


2.2 


352.2 


5-10 % 


4.04 


294.7 


10-15 % 


5.27 


245.6 


15-20 % 


6.26 


204.2 


20-30 % 


7.48 


154.5 


30-40 % 


8.81 


103.8 


40-50 % 


10.01 


64.9 


50-60 % 


11.11 


36.6 


60-70 % 


12.06 


18.8 


70-80 % 


12.96 


7.5 


80-90 % 


13.85 


4.4 



its energy as described in the previous sections. During the time of its passage strings partially 
decay and loose their colour. This effect is taken into account according to Eq. ([9]). The average 
length l n k traveled by the particle emitted from the k-th cluster of n strings depends both on the 
distribution of clusters and on the direction of the emission. Due to azimuthal asymmetry of the 
cluster distribution following from the asymmetry of the overlap area the average distribution of 
emitted particles will depend on the azimuthal angle and lead to non- vanishing anisotropic flows. 

In the simulation of an event one has, first, to determine the distribution of clusters in the 
overlap area and, second, for each cluster to find the average length of the path which the 
emitted particle has to travel inside the cluster matter. The final distribution of emitted particles 
in the transverse momentum and azimuthal angle will be given as a sum 

P(p,<P) = Cj2^nke (12) 

n,k 

and the distribution in the azimuthal angle only 

P{4>) = C j ' dp 2 P{p^). (13) 

(We use the thermal distribution ([3]) to be able to move into the region of p of the order of several 
GeV/c) 

With the number of strings N not very large (below 100) the described procedure is realizable 
on a computer for a reasonable processing time. However this time grows very fast with the 
number of strings (as approximately ~ N 3 ). This motivated a somewhat simplified approach to 
our simulations. Instead of taking clusters at different locations in the transverse space and of 
different geometric forms, we assumed them to form a square lattice in the transverse space with 

the side length a = ^nr^. Throwing of discs transforms these primitive structures into clusters of 
strings (" ministrings" ) with a variable number of fused simple strings equal to the number of disc 
centers found in a given lattice cell and thus with a variable cell-dependent percolation parameter 
p. In this manner we avoid the study of all complicated geometric structures which arise when 
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a large number of simple strings overlap. Calculations of the momentum distribution of emitted 
particles in such a picture shows that it models the actual cluster formation quite successfully [37] . 

Note that in Eq. (TTZj) it is assumed that multiplicity fi n k of the cluster is fixed by its area and 
colour according to Eq. (|10p. This multiplicity is the average over events. On the event-by-event 
basis one may take the multiplicity distributed according to Poisson's law around this average, 
as advocated in [38J, where it was shown that such distribution leads to the negative binomial 
distribution of observed particles. We shall see that such stochastic emission does not change 
significantly the final results for the flow coefficients although naturally somewhat complicates 
numerical calculations. 

Note that the fluctuations in the distribution of strings and multiplicities are not the only ones 
in the dynamics of particle production. Also the numbers of participant nucleons N par t and strings 
inside each of them N s t r fluctuate as well as the values of the impact parameter b. However these 
additional fluctuations seem to have little to do with the azimuthal anisotropy, as supported by 
our results, which show that fluctuations in the multiplicity do not influence the final anisotropy. 
Therefore in our calculations we have not taken into account fluctuations of N par t, N s t r nor b, 
although they can be essential in other context [39] 

In our model, with the coefficient a for the string color decay fixed at the value 0.03, as explained 
in the previous section, only a single dimensionless parameter k remains, which characterizes the 
loss of energy in passing through the string field and is to be extracted from the experimental 
data. It was adjusted to fit the experimental value of Vi integrated over the transverse momenta 
for Au-Au collisions at 200 GeV in mid-central events. 

On the event-by-event basis the final flow coefficients v n are found from the distribution of 
emitted particles P(p, <fi) at a given p or P(<fi), integrated over all p by the standard formulas. Let 
in a given event 

a n = J d(j>cos(n4>)P(4>), b n = J d4>sm(ntp)P(<f>) (14) 

then 



v n = V< + bl ,n = l,2,.... (15) 
a 

4 Numerical results 

Our numerical calculations were performed for Au-Au collisions with the profile function taken ac- 
cording to the Saxon- Woods formula and cluster multiplicities distributed according to the Poisson 
law. The ordinary string tension parameter To was taken to agree with the slope of the spectra in 
the soft region: To = 0.08 (GeV/c) 2 The adjusted value of the quenching coefficient turned out to 
be k = 0.058. The number of simulation was taken 100. Our Monte-Carlo results obtained for v<i 
and integrated over the transverse momenta and averaged over events are shown as a function 
of b in Figs. [2]|4]for Au-Au collisions at energies 62.4, 200 and 2760 GeV respectively. In these and 
the following figures error bars show fluctuations of the flow coefficients around their average values 
from event to event, not the insufficient precision of the calculations due to the finite number of 
simulations. One observes that these event-by-event fluctuations are quite large. 

To show higher harmonics we present v n for n = 1, ...8 integrated over pt at centrality 40—50 
% in Figs. [SJ- [7J again with their event-by-event fluctuations from the average. 

Finally in Figs. [8l-llllwe present our results for the pt dependence of V2 and 1*3 at two centralities 
10-15 % and 40-50 % for Au+Au collisions at energies 200 GeV and 2.76 TeV 
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Figure 2: V2 and i>3 integrated over pr as a function of b for Au-Au collisions at E = 62.4 GeV. 
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Figure 3: Same as for Fig. E3 at E = 200 GeV 
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Figure 4: Same as for Fig. [2] at E = 2.76 TeV. 
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Figure 5: v n , n = 1, ...8 integrated over px at centrality 40—50 % for Au+Au collisions at 62.4 GeV 
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Au+Au at 200 GeV 
Centrality 40 - 50 % 



Figure 6: Same as for Fig. [5] at E = 200 GeV 



Au+Au at 2.76 TeV 
Centrality 40 - 50 % 



Figure 7: Same as for Fig. Oat E = 2.76 TeV 
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Au*Auat 200GeV 



Centrality 30 - 40% 



Am A j at 2.76 TeV 



Centraity 10-15% 



Figure 8: V2 and i>3 as a function of px for Au+Au collisions at 200 GeV and centrality 10—15 % 



Au*Au at 2.76 TeV 



Au+Aj al 2.76 TeV 



Figure 9: V2 and v% as a function of pr for Au+Au collisions at 200 GeV and centrality 40—50°-% 



Centrality 30-40% 



Figure 10: V2 and v% as a function of pr for Au+Au collisions at 2.76 TeV and centrality 10—15 % 



Figure 11: V2 and v% as a function of pr for Au+Au collisions at 2.76 TeV and centrality 40—50 % 
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Inspecting these results we see that the b- and pr dependence of v 2 as well as their absolute 
values well agree with the experimental data (which of course is based on the proper adjustment of 
coefficient kq at E = 200 GeV and medium centrality). The behaviour of the triangular coefficient 
t>3 is also found to agree with the experimental observations. Its average values are calculated to 
be somewhat smaller than the experimental values. However their event-by-event fluctuations from 
the average turn out to be very large, so that taking them into account may put the experimental 
values well within the calculated ones. 

To analyze the important components of fluctuations leading to relatively large values of V3 we 
repeated our calculations substituting Poisson distributed cluster multiplicities by their average. For 
Au+Au collisions at 200 GeV and centrality 40—50 % this gives values for v 2 and v% integrated over 
p T 0.594E-01 ± 0.145E-01 and 0.574E-02± 0.403E-02 respectively to be compared with 0.588E-01 
± 0.186E-01 and 0.668E-02 ± 0.399E-02 with Poisson distributed multiplicities. As one observes 
the difference is insignificant. On the other hand the importance of quenching is found to be 
overwhelming. Taking kq = but keeping the Poisson distributed multiplicities one obtains values 
for v 2 and v 3 0.180E-02± 0.839E-03 and 0.112E-02± 0.622E-03, that is practically zero. 

5 Conclusions 

We have performed detailed Monte-Carlo simulations to study anisotropic flows in the percolating 
string scenario. We have confirmed that the colour string model with fusion and percolation can 
successfully describe the observed elliptic and triangular flows in high-energy heavy-ion collisions. 
An important ingredient in this description is anisotropy of the string emission spectra in the 
azimuthal direction which follows from quenching of the emitted partons in the strong colour field 
inside the string. As a mechanism for this quenching we used the radiation energy loss during 
propagation of a fast charged particle in a constant field, borrowed from the QED. Upon adjusting 
the parameter of quenching, this allowed to describe the data quite well both in their energy 
centrality and transverse momentum dependence. We have also studied higher anisotropic flow 
coefficients from to vs which are found to be small as compared to v 2 . 

To compare with our earlier calculations in the grossly oversimplified picture in [13| we have 
found more pronounced energy dependence. This follows from the explicit dependence of the 
quenching in the QED on the field strength, which is translated into dependence on the string 
tension in the colour string picture. So, in contrast to our phenomenological formula adopted in 
|13j . this quenching is not purely geometrical but grows with the percolation parameter even when 
all strings are fused into a single cluster occupying the whole overlap area. 

Our results are based on strings which have infinite dimensions in rapidity. So from the start 
they refer to the central region of nearly zero rapidity and thus do not allow to study the rapidity 
dependence of the flow coefficients. To do this we have to introduce strings of finite rapidity length 
and in this way take into account energy conservation. This complicates our picture substantially 
and will be the object of our further studies. 
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